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PREDICTION  OF  LIQUID -YAPOR  EQUILIBRIA  FROM  AN  EMPIRICAL 
EQUATION  OF  STATE  FOR  MIXTURES* 

by 

George  Galloway  Dauplaise  Town 
I.  INTRODUCTION 

This  report  set*  forth,  a method  of  calculating  the  equilibrium  phase 
boundary  for  which  one  needs  only  an  equation  of  state  for  mixtures  and  the  set 
of  accompanying  constants  for  each  of  the  pure  components. 

The  advent  of  high  speed  automatic  calculating  machines  has  made 
feasible  the  extensive  calculations  necessary  to  determine  the  pressure, 
volume,  temperature  and  mole  fraction  properties  of  a multicomponent  system 
at  the  liquid -gaseous  phase  boundary.  Connolly,  Frankel,  and  Sage*  have 
developed  a method  of  determining  the  phase  boundary  utilizing  the  International 
Business  Machine  Electronic  Calculating- Punch  604.  This  method,  however, 
uses  experimental  data  a©  initial  guesses,  making  a calculation,  and  then 
proceeding  in  an  iterative  process.  The  data  is  available  for  only  a relatively 
few  of  the  possible  systems  and  therefore  the  method  is  of  limited  applicability. 


T.  J.  Connolly,  S.  P.  Frankel,  and  B.  H.  Sage,  Report  No.  3036, 
American  Documentation  Institute,  Washington.  D.  C.,  (1951). 

* This  work  was  carried  out  at  the  University  of  Wisconsin  Naval  Research 
Laboratory  under  Contract  N7onr-2S5ii  with  the  Office  of  Navo.1  Research. 
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II.  THEORY 


It  has  beer  shown  by  the  arguments  of  thermodynamics"  that  the 


necessary  and  sufficient  conditions  for  equilibrium  are: 


PG  ' Pl. 
TG  * TG 


Where: 


A iG  = iJ 


(i  = 1,  2,  ..  . n ) 


Y* I 

■ i 


P is  the  pressure  cf  the  phase. 

T is  the  temperature  of  the  phase. 

/L  is  the  chemical  potential  of  the  ith  component. 

The  subscripts  G and  L refer  to  the  gaseous 
and  liquid  phases  respectively. 

From  this  set  of  equations  the  P-V  T and  mole  fraction  properties  of  the 
liquid  and  gaseous  phases  at  the  boundary  can  be  uniquely  determined  if  an 
equation  of  state  for  mixtures  is  assumed,  the  chemical  potential  being 
evaluated  from  the  equation  of  state' . 

The  particular  equation  of  state  used  was  that  devised  by  Benedict, 

4 

Webb,  and  Rubin  . The  form  of  the  equation  given  below  refers  to  one  mole 


J.  W.  Gibbs,  ''Equilibrium  cf  Heterogeneous  Substances",  Collected  Works, 
5b,  Longman,  Green  and  Co.,  New  York,  (1928). 

J.  A Becitie,  Phys.  Rev.,  36,  132  (1930). 

M.  Benedict,  G.  B.  Webb,  and  I,  C.  Rubin,  J.  Chem.  Phys.,  8,  334, 
(1940);  10,  741  (1942). 
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**£  substance: 

p 


XTd  4-  {BqST  - Ao  - Co/T2)d2  + (b»T  - a)** 
4-  a*  d6  i (c/T2)d3  1/1  4 1 d2)  e“yd2_/ 


.3 


U> 


in  which: 


P is  the  pressure  in  atmospheres, 

K.  is  the  gas  constant, 

T is  the  temperature  in  degrees  Kelvin, 
d is  the  density  in  mole  per  liter. 

The  constants  3>0  , , h , a , c , , 2r  _ are  determined  from 

the  constants  for  the  pure  components  - BQi  , AQi  , , aj  , , 

- an.d  their  respective  mole  fractions  (xj)  in  the  following  way*: 
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The  chemical  potential  determined  u*ing  equation  2 i»: 

A - m°  -RTInfdXi  RT)  + df'4RT  ExjfCBo'j ’*  +(Bo V,14]* 

-2[4Mji!w-rvT2)[e0fc^'4} 
+ i dJ  [ RT  (b*b|Vte  - (&*  a\V*] 

Tid![a(c(!  *,y'~  + fT  *,V/S'j 

+ 3(C2Ci)!ft  Ai  ri-e""12  e'rai 

T*  L 5 

_ ML  cuf  A.  1-e"*1  - e -r*'  -a 


VVFJ 


in  which  J4  is  the  chemical  potential  of  the  ith  component  at  PQ  and  Tq  of 
the  standard  state. 

M ± * defined  below,  was  used  in  the  calculations  and  givesthe 
deviations  from  ideality  of  the  particular  system; 


A -J\  -V-RT!hfRT) 


IE,  SYSTEM  STUDIED  AND  METHOD 


The  binary  system  of  n-beptane  (component  1)  and  nitrogen  (component  2) 
was  investigated.  The  sets  of  constants  used  ior  the  pure  substances  are  listed 
in  table  1.  The  mole  fraction  of  n-  heptane  is  designated  as  Ks.*1  and  the  mole 
fraction  of  nitrogen  as  "i  - x”  since: 


£x.  =1  (i.  = 1,  ?,) 
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TAB  JOE  I 

The  r«nstant6  for  nitrogen  and  n-Hept&nc  used  in  the  Benedict- 
Webb-Rubin  equation  of  a i-ale. 


5 

n~ heptane 

0. 199005 

1 1 

X k • w wn  V 

4. 7457' 

U . i 5 i 9 54 

: <*!  ~ V 

•7  4 /%  o 

0.0043561 1 


0.0458000 
) . 19250 

C . 05S8907 

0 . u‘U  i9c  I 94 
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° M.  Benedict,  G.  B.  Webb,  JL.  C.  RnMn„  <m«  L-.  t?'ri»-ad,  Chem 
Eng.  Prog.,  47,  421,  table  1 (1951). 

A private  communication  from  D r.  bianaou  Benedict.. 
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To  solve  the  set  of  equations: 
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T = T 
G L 


G ~ L 

/iG  = 

/ 


(7) 

<8) 


i 9) 


iL 

>*2G  = ^2L  (i0> 

T and  P were  chosen  equal  across  the  phase  boundary.  The  pair  of  mole 

fraction  of  n-hepiane  (x^-.  and  x^)  which  satisfy  equation  (9)  and  (10)  at 
the  given  T and  P were  determined  by  the  method  of  repeated  plotting  on  a 
larger  scale ^ of  and  w'  (of  both  phases)  as  functions  of  and 

Sets  of  x's  were  chosen  fez  both  the  Liqujri  and  vapor  phases  which  bracketed 
the  answers.  These  x‘s  were  then  used  to  compute  and  ^or  the 

P and  T under  consideration.  The  approximate  solution  of  the  simultaneous 
equations  was  determined  from  the  graph  by  inspection.  These  values  of  x^ 
and  x^  were  then  used  to  choyii.  seta  of  more  closely  spaced  x’s  which  were 
in  turn  ut-ed  to  evaluate  the  and  /<'  at  the  Z'Vr:''  ? and  T.  Figure  ] 

shows  me  last  approximation  to  and  at  49d.  i30  and  IOC  alms. 


j i ' . y 

j and  /* were 


aicuiated  with  the  o»  the  in  ;er  national  Business 


Machine  Model  1 Card-programmpH  Electiomc  Calculator  for  a given  x , P and 
T frem  equations  (2)  , (3)  , and  (4)  The  programming  deck  of  cards  used  in 
the  calculator  first  deti*rmif>?d  th?  density  (d)  from  equations  (2)  and  (3)  by 
four  iterations  oi  the  Newton  Raphcnn  ’^etasi"'.  This  value  of  d , along  with 


equations 


the  giver,  x and  T , was  then  used  to  determine  and  J*' 

(3)  and  (4). 

7 

J.  B Scarborough,  Numerical  Mathematical  eralyais,  p.  190,  The  Johns 
Hopkins  Press,  Baltimore  (1950). 


Ibid.,  p.  20  3. 
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IV.  RESULTS  AND  DISCUSSION 


Four  temperatures  were  investigated:  298.  13°K,  373,  13°K,  448.  13°K, 

498.  13°K.  The  results  are  plotted  in  figure  2 and  bated  in  table  2.  The 

results  at  29S.  13°K  were  compared  with  the  experimental  results  of  Palmer, 

,a  Q 

HirecMeider,  and  Boyd  at  295.6  K'.  The  comparisons,  listed  in  table  4, 
were  made  by  plotting  In  (x./  as  a function  of  1/T  for  a particular  pressure 
and  extrapolating  to  1/295.  6°K.  The  probable  error  of  :che  calculated 
compared  to  the  experimental  results  is  0.004.  The  experimental  r«  suits  of 
Boomer,  Johnson  and  Piercey  for  the  n-heptane  nitrogen  system  are  in 


disagreement  with  the  results  of  this  work.  The  mole  fraction  of  n heptane 
in  the  vapor  phase  at  298.  13°K  obtained  by  Boomer,  Johnson  and  Piercey10  was 
0 0104  as  compared  to  G.  002160  for  the  calculated  result. 

The  range  of  reliability  of  the  calculated  x^-.  at  298.  13°K  is  from  30  to 
250  atm.  This  limitation  is  caused  by  equation  (2)  failing  at  densities  greater 
than  twice  the  critical  density.  This  failure  is  the  reason  that  the  upper  portion 
of  the  log  P versus  log  x carve  could  not  be  closed  and  the  gac  portion  of  this 
curve  could  not  be  extrapolated  to  the  experimental  vapor  pressures11  of  n-heptane 


q 

In  j^.  rainier,  J.  O.  H.rschfeider,  and  C.  A.  Boyd,  Univerehv  of  Wiaco«ci 
rla . u..  HciucZvit  i^ano rater y,  c li"  * « « , 1952.). 

[ f\ 

E.  H.  Boomer,  C.  A Johnson  and  A.  G.  A.  Piercey,  Can.  J.  Res  , 

16B,  396  (1938). 

S.  Young,  J.  Chem  Soc.,  73,  675  (1898). 
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TABLE  2 


The  results  of  the  calculations  for  the  phase  boundary  of  nitrogen  and  n-heptane 


M\  .. 


P (atm)  XG  x 100 

x 100 

dG 

298.  13° 

dL 

K 

( 1 - atm /mole)  ! 

1 1 - atm/  xno 

0,2 

21.63 

99.974 

0.008186404 

6.68404 

-155. 313 

-123. 70 

0.5 

8.70  3 

99.922 

0.020456330 

6.686110 

-155. 294 

- 97.38 

1 

4.  392 

9". 844 

0.040910 

6.6893 

-155,  23 

- 79.33 

5 

0.9414 

99.  19 

A i efA 

U . L\J  'i  l U 

b.  7 16 

-154. 79 

■7  n i * 
- j 7 . x x 

10 

0. 5131 

98.  38 

0.40981 

6.749 

-154.  24 

- 22. 07 

30 

0. 2402 

95.  16 

1.23272 

6.883 

-152.02 

4.  79 

50 

0 . 200  2 

91.98 

2.05560 

7.025 

-149. 79 

17.  22 

70 

0. 1966 

88.85 

2.87149 

7.  145 

-147. 53 

25.44 

100 

0, 2160 

84.  20 

4.07868 

?„  383 

-144.  12 

34.  1 6 

150 

0.  2862 

76.  56 

5.99793 

7.777 

-138, 31 

44.  20 

170 

0. 3246 

73.  52 

6. 72485 

7.945 

-136. 35 

46.97 

20  Q 

0.  3923 

69.00 

7.76493 

8.  209 

-132.43 

51.49 

250 

0. 5306 

61.60 

9. 35527 

8.682 

-126. 45 

57.  34 

300 

0. 6973 

54.45 

10.76187 

9.  197 

-120. 39 

62.  31 

400 

1.090 

42.  45 

13.04954 

10. 250 

-107.96 

70. 653 

500 

1.489 

34.85 

14. 7440 

11. 149 

- 95.01 

77. 760 

600 

1.810 

30. 88 

16.0136 

11.806 

- 3 1 . 50 

84  13 
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TABLE  2 (continued) 


P (atm) 

XG  x 100 

Xh  x 100 

G 

dL 

( 1 -aim/mole) 

>2 

( 1 - atm/ mole) 

373 

I j XV 

!• ! 

! 1 

•y 

U 

55,  10 

99.812 

0.0675290 

5.93204 

-104  . 30  3 

- 107.  30 

3 

37.43 

99.61 

0.  100408 

5.93887 

-104.  197 

- 79.9 

• 

S 

J 

5 

23.06 

99.  21 

0 . 16626 

5.9525 

-103.,  98 

- 63.  1 

j 

5 

i 

10 

12.  20 

98.  20 

G.  330  29 

5.9869 

-103.45 

- 38.00 

* 

* 

» 

• | 

25 

■“  / 

D.  O IV 

r i r»f* 

7">.  A03 

0. 8390 

6.  119 

-10  1.  84 

7,  74 

* | 

50 

3.612 

90.18 

1,63033 

6. 2740 

- 99. 142 

14.  31 

70 

3.  113 

86.  17 

2. 2651 

6.427 

- 96.98 

24.  94 

i 

100 

2.864 

80.  14 

3.  2043 

6, 6760 

- 93,73 

36.  24 

1 

it 

150 

2.957 

69.88 

4.699 

7.  117 

- 88. 33 

49.  24 

| 

200 

3.  355 

59.  25 

6.  1001 

7.634 

- 82.95 

58.  67 

k 

j 

250 

4.021 

48.  19 

7. 3980 

8.  232 

- 77.64 

66.  17 

1 

i 

. 

300 

5.003 

37.  11 

8.  584 

8.907 

- 72.410 

72.478 

ICO 

6.  568 

26.  82 

9.6445 

9.635 

- 67.215 

78.0073 

400 

10.  5 

16.0 

10. 526 

10.451 

- 61.96 

82.90 

4 
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TA.&LE  2 (continued) 


P (atm) 

"VT  _ y 100 
'u 

X.  x 100 

dG 

dE 

A 

( 1 aim/ mole) 

/ 

A 

( 1 - atm /mole) 

448 , 

1 3°K 

8 

82.  69 

99.  59 

0.25804 

5. 16577 

- 70.850 

-144.9 

10 

68.90 

99. 11 

0. 31428 

5.  17830 

- 70.  643 

-102.  8 

12 

59.  31 

98.61 

0. 37042 

5.  19094 

••  70.435 

- 70.05 

15 

49.39 

97.88 

0.45459 

5.  210 

- 70.  120 

- 54.40 

25 

32.90 

95.435 

0.7343 

5. 2740 

- 69.08 

- 26. 15 

cn 

M*V 

20. 28 

89.  30 

1.4272 

5.4382 

- 66.49 

5.  35 

Kl/H 

IU 

16.95 

84.  37 

1.9745 

5.  574 

- 64.40 

19.42 

100 

14.98 

76.78 

2.  7827 

5.486 

- 61.31 

3j.  7 ( 

150 

14.93 

63.  20 

4.099 

6.  153 

- 56.23 

49.61 

200 

18.  13 

46.  19 

5.415 

6.480 

• 51.31 

62. 00 

220 

22.  21 

35.  56 

6.014 

6.  506 

- 49.404 

65.905 

498. 

13°K 

20 

82.  28 

98.420 

0. 68650 

4.  4991 

52. 830 

- 86.55 

50 

47.  44 

89.  20 

1. 509 

4.  625 

- 49.80 

- 8.63 

100 

39.98 

70.  69 

2.  999 

- /.AC 

• u / 

A A O 1 

••  ~*I . ox 

0 An 

£ / . 77 

110 

42.08 

64.91 

3.  3S45 

i.  60  6 

- 43.8314 

32.643 

115 

44.66 

60.  28 

3. 6456 

4.4687 

- 43. 3573 

34. 7737 
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TABLE  3 


Vapor  Pressure 

of  n-heptane 

T°K 

Pressure  (atm.) 

298.  13 

0.0612 

373. 13 

1.046 

448. 13 

6.099 

498.  13 

14.  38 

i i 


i 


i 


* 

. r 

TABLE  4 

Comparison  of  calculated  with  the  experimental  results 

of  Palmer,  Hirschfelder , and  Boyd  at  295.  6°K. 


(atm. ) 

Cal.  r 100 

(j 

Exp.  xG  x 

150 

0.258 

0. 255 

100 

0.  194 

O IOC 

V • * / »> 

70 

0,175 

0.  130 

50 

0.  179 

0.  184 

30 

0.217 

0.  227 

1* 

ii 


i'  . • 
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as  the  other  cui  ves  do. 

I he  graphical  method  o£  solution  developed  could  be  used  to  obtain 
any  degree  of  precision  desired,  but  this  would  not  be  practical.  Rather 
the  results  of  the  graphical  method  would  be  used  as  ihc  initial  guesses  in  an 
iterative  method  such  as  the  Nev/ton-Raphson  method  for  simultaneous  equations. 
However,  due  to  the  limited  accuracy  of  the  Benedict- Webb -Rubin  equation  of 
state  such  extensive  calculations  are  not  warranted. 


V.  CONCLUSIONS 

The  following  conclusions  can  be  drawn  Irom  this  work: 

1.  A practical  method  has  been  developed  by  which  one 
can  determine  the  pressure-mole  fraction  curve  of 
the  liquid-gaseous  phase  boundary  for  a given  temperature 
f'om  an  empirical  equation  of  state. 

Z.  The  calculated  result*  from  our  equation  of  state, 

using  the  constants  of  the  pure  components  in  the  equation, 
are  in  agreement  with  the  exper  imental  results  of  Palmer, 
Hirschfelder,  and  Boyd,  but  are  in  disagreement  with  the 
experimental  results  of  Boomer,  Johnson  and  Piercey. 


- e ciisi.  «. 

4>* 


‘■t'S* 


